Jarzynski equation for a simple quantum system: Comparing two definitions of work 
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The validity of the Jarzynski equation for a very simple, exactly solvable quantum system is 
analyzed. The implications of two different definitions of work proposed in the literature are inves- 
tigated. The first one derives from measurements of the system energy at the beginning and at the 
end of the process under consideration making the work a classical stochastic variable with transi- 
tion probabilities derived from quantum mechanics. In the second definition an operator of work is 
introduced and the average in the Jarzynski equation is a quantum expectation value. For the first 
definition a general quantum mechanical version of the Jarzynski equation is known to hold. For 
the second one the Jarzynski equation fails to yield the free energy difference at low temperature. 
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INTRODUCTION 

The fluctuation properties of physical systems driven 
strongly away from thermal equilibrium are of central 
importance in statistical mechanics. Recent progress in 
this field gathers from a handful of exact and remark- 
ably general theorems characterizing the properties of 
exponentially rare trajectories. Although closely related 
to each other it became customary to divide them into 
two groups: non-equilibrium work theorems relating the 
work done in an irreversible process to the free-energy 
difference between initial and final state [H, B S 0] j an d 
non-equilibrium fluctuation theorems concerning the ra- 
tio between the probabilities of entropy-producing and 
entropy-consuming processes [a, |fj|, |7|, |8|, |fj, [Kj, . 

Several proofs of these theorems are by now avail- 
able if the microscopic dynamics is described in terms of 
classical mechanics. Particular examples include micro- 
canonical Hamiltonian dynamics d EH , Hamiltonian d\ 
namics coupled to different types of thermostats 
and stochastic dynamics of Langevin type [7|, Ifj, [l( 

The situation is less satisfactory for the case in which 
the system is described by quan tum mechanics at the 
microscopic level [II El 0, El El El El El III IE , 
I22I [23| . Central problems in generalizing non-equlibrium 
theorems of the described type to the quantum realm 
are the definition of the work performed on the system 
in the irreversible transition and the identification of the 
proper average to be taken. If the average is interpreted 
as an average over quantum mechanical transition proba- 
bilities for the evolution between initial and final Hamil- 
tonian HI EI El El E3| the description of the quantum 
mechanical system is formally similar to a classical one 
and quantum versions of the fluctuation theorems can be 
proven in analogy to the classical case. Alternatively the 
work may be related to the difference of the Hamiltonian 
of the system in the Heisenberg representation at differ- 
ent times [lE3,0,[2l- The average to be performed is 
then a quantum mechanical expectation value with the 
initial density operator. 

In the present note we investigate the validity of a lim- 



iting case of the Jarzynski equation for a very simple 
quantum system. This system is exactly solvable and 
physically transparent and the implications of both defi- 
nitions of work and corresponding averages can be stud- 
ied analytically. We will show that both approaches are 
compatible with the second law of thermodynamics and 
correctly describe the classical limit f3h — > but yield 
different results in the quantum regime. We will indicate 
the formal as well as the physical reason for the failure 
of the Jarzynski equation in the second approach which 
we nevertheless believe to be the more appropriate one. 



THE MODEL 

We consider a one-dimensional quantum harmonic os- 
cillator of mass m and use units in which m = h = 1. 
The Hamiltonian is given by 
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J ~2 
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with the usual meaning of the symbols. The system is 
supposed to be at equilibrium with a heat bath of inverse 
temperature /3. From elementary quantum statistics the 
free energy is given by 
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The system is now decoupled from the heat bath and 
the frequency of the oscillator is changed from u>q to uji 
according to 



t 



if 
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where tf denotes the final time of the process. For < 
t < tf the system is hence characterized by a unitary 
time evolution described by 



U(t) = T> cxp ( - - / dt'H(t') 



(4) 
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induced by the time-dependent Hamiltonian 



H(t) 
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where T> denotes the time ordered product. The classi- 
cal version of this system was discussed in [25j, |26j . 

Since there is no coupling to the heat bath during the 
process the work W necessary to accomplish the change 
luq uji equals the change in energy AE of the system. 
The Jarzynski equation [l| hence stipulates 



-0AE } = e - 



f3AF 
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where (...) denotes the average over repeated realizations 
of the process and 
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as follows from 



MEASURING ENERGY TWICE 



The definition of work advocated in 0, 0, [H, El | 
derives from a measurement of the energy of the system 
at the beginning and at the end of the process. Mea- 
suring H at t = yields an eigenvalue E^ of Hq with 
probability 



p(0) _ e /3(F -£< 0) ) 



and leaves the system in the corresponding eigenstate 
\ipn ) of Hq. At time tt the system is hence in state 
U (t f)\ipn^ ) and measuring H anew yields an eigenvalue 
Em" 1 of Hi = H(tf) with probability 
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whith | 1pm ) denoting the corresponding eigenstate of Hi. 
The energy difference AE is now identified with Em — 
E^ and the average of a function f{AE) is defined as 

(f(AE)) = Pr^ E W ™ - E&) (10) 

n m 

One then easily finds 0, Q (H 

{e -^E ) = J2p(P)J2W mn e-^ E (11 

n m 

= e ^ (0) ^e-^ 0, |(^ ) |t>(i/)l^ 0) )| 2 e-^^ , -^ ) ) 
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Hence for this definition of work the Jarzynski equation 
can easily be extended to quantum systems. The second 
law in the form (AE) > AF then follows in the usual 
way from Jensen's inequality. 

In this approach no genuine quantum observable of 
work associated with a self-adjoint operator is intro- 
duced. Instead a classical energy difference is combined 
with quantum-mechanical transition probabilities. This 
makes the proof of the Jarzynski equation very similar 
to the case of a classical stochastic process described by 
a Master equation [25| . In fact all that is needed from 
quantum mechanics is J2 n W mn = 1 following from the 
unitary time evolution 



OPERATOR OF WORK 

In the present setup the work performed is equal to 
the energy difference AE of the system. Since the (time- 
dependent) energy of a quantum mechanical system is 
associated with the Hamilton operator in the Heisenberg 
representation it is reasonable to represent AE by the 
Hermitian operator @, H 0, 01 



AE = U\t f )H(t f )U(t f )-Ho 



(13) 



The relevant average of a function f(AE) is then given 
by 

{f(AE)) = Tr (p f(W(t f )H(t f )U(t f ) - H j) , (14) 

(g) where po denotes the equilibrium density matrix corre- 
sponding to H and Tr stands for the trace. Eq. (fT4"|) rep- 
resents a genuine quantum average over the initial con- 
dition of the process. It characterizes the time evolution 
of the energy of the system between t — and t = tt 
and does not include contributions resulting from mea- 
surements. 

It is easily shown that the two averages defined in (fT0|) 
and (fT4"|) respectively give identical results for functions 
f(x) involving only linear and quadratic terms. There- 
fore both procedures yield the same result for (AE) and 
are hence both consistent with the second law. However, 
whenever f(x) contains higher order non-linearities, as 
for f[x) — e~ x , the results are only the same if the 
Hamiltonians corresponding to different times t and t' 
commute with each other, [H (t), H(t')] = 0. In this case 
we get indeed 



P pm Tr ( p -PH a p -f3(Ui(tf)HiU(t f )) RH a 



-pAF 



and the Jarzynski equation again holds. However, in gen- 
eral [H(t),H(t')} ^ 0, therefore 
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and the Jarzynski equation will not be fulfilled. 

We note that the proof of a quantum analogue of the 
Jarzynski equation given in (l2| starts out with a defini- 
tion of a work operator similar to (|13|) . In the subsequent 
proof, however, an operator of exponentiated infinitesi- 
mal work is introduced which is somewhat ambigous. It is 
argued that this ambiguity becomes irrelevant in the limit 
tf — > oo. For general tf, and in particular in the limit 
tf — > we are interested in here (see below), however, 
this ambiguity disappears only if [H(t), H(t + dt)] = 0. 



EXACTLY SOLVABLE LIMIT CASE 

For the toy model introduced above the implications 
of the second definition of work can be analyzed exactly 
in the limit tf — > 0. The corresponding limit case of 
the Jarzynski eq uation is also known as thermodynamic 
perturbation [281 ] . In this case of an instantaneous switch 
from Hq to Hi we have 



U(t f ) = + H 1 )t f + 0(t 2 ) 
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and therefore 



AE = Hi - H 
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The work operator is hence diagonal in the position re- 
spresentation and using for the density operator [2^ | 



p (x,x ) = y — 



tann 



exp ^ — ^ ((x + x') 2 tanh — (x — x') 2 coth 

(18) 

a simple calculation gives for the Jarzynski estimate A J 
of the free energy difference 
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The average of AE itself is given by 



(AE) 



0UJ O 



^^°coth — 
4uj 2 



(19) 



(20) 



Comparison of ([19]) and (|20p shows that for the present 
system we have the exact relation 
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From and H| we infer (AE) > AF for all with 
the equality holding only for u>i — luq. This is in ac- 
cordance with the second law of thermodynamics and 
the fact that the instantaneous change of w is a non- 
equilibrium process. The inequality holds down to T = 
where we have 



(AE) = (^ 
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with the inequality sign resulting from the Rayleigh-Ritz 
variational principle. 

On the other hand the results |0 and (|T9|) differ from 
each other for all values of 0. Therefore the Jarzynski 
equation does not hold for the system under considera- 
tion if the second definition of work is implemented. In 
view of the discussion in the previous section this is in 
accordance with [Hi, Ho] ^ 0. 

In the high-temperature limit, — > 0, we find 
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such that the leading behaviour is the same in accordance 
with the validity of the Jarzynski equation for classical 
systems. The first differences occur already in the O(0 2 ) 
corrections. Since the expansion in is equivalent to an 
expansion in h we find that in the classical limit devia- 
tions from the Jarzynski equation occur at order 0(h 2 ). 
A similar result was obtained recently for quantum cor- 
rections to the transient fluctuation theorem [2lj |. 
In the opposite limit, — > oo, we have 
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showing qualitative differences between the two quanti- 
ties in the quantum regime. Whereas AF tends, as it 
should, to the difference between the respective ground- 
state energies of Hi and Ho the estimator A J approaches 
zero (as long as u>i > ljq). The situation is illustrated in 
Fig.ffl 

The differences between the two definitions of work 
for the system under consideration become particularly 
clear when considering the case u>i < u>o, as illustrated 
by the inset of Fig. [TJ Indeed, from (fj"9|) we infer that 
A J diverges at = C where C is the solution of the 
equation 



1 , ,/&<*>(), 
^ tanh( ^- } 



2ujo 



(26) 



The origin of this divergence is easily understood intu- 
itively. For u>i < loo the system gains energy by the 
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FIG. 1: Free energy difference AF as given by (J?) (full line), 
Jarzynski estimate A J defined in (|19p (dashed), and average 
work {AH)o given by (|20|l (dotted) as function of the tem- 
perature T = 1/(3 for o>o = 1 and ui\ = 2. The inset shows 
the same quantities for uo — 1 and wi = 0.5. 



change in u> and accordingly AE as denned in (fl~7|) is neg- 
ative semi-definite. The system receives a large amount 
of energy if the particle is far from the origin when the 
change in lo occurs and these contributions may induce a 
divergence of (e~@ AE ) for sufficiently large 0. In the clas- 
sical case no such divergence may occur since the initial 
canonical distribution po(x) ~ exp(— 0uJqX 2 /2) confines 
the particle with increasing to a, smaller and smaller 
vicinity of the origin which exactly counterbalances the 
impending divergence. In the quantum case and using 
the second definition of work, however, we find from (I18|) 
in the limit — > oo 



work it is worthwhile to ask whether there is an operator 
A such that 



(e-P A ) = e-? AF . 
This operator has to fulfill 

T[e -pA e -fi6 0=T[e -fi6 1 

and one possible choice is 0, HI] 
1 



(28) 



(29) 
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Explicitly we find in the present case 
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where the real coefficients a, b and c can be determined 
from the Baker-Hausdorff-Campbell [30j formula. The 
first terms are of the form 
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The operator A given by (f3"Tj) does satisfy ([2^)1 but does 
not represent a physical observable since it is not Hermi- 
tian. Moreover it is rather specific for the present exam- 
ple. 



pa[x,x) J — 



l^ 0) (-)l 2 



(27) 



and the size of typical excursions of the particle from 
the origin becomes independent of 0. Correspondingly, 
for sufficiently large values of the Jarzynski average 
{e~ l3AE ) will diverge. It is also instructive to consider 
why this divergence does not occur in (|11| . The rea- 
son the sum over n cannot diverge is the same as in 
the classical case: large values of En"* are suppressed by 
correspondingly small values of Pn . The m-sum can- 
not diverge either since for given n the energy difference 
AE = Em^ — E^ becomes positive for sufficiently large 
m in spite of u)\ < Wq. In other words, even though the 
energy difference is negative semi-definite the measure- 
ment of the energy at the end of the process may yield 
higher values than the one at the beginning which renders 
the expression (jlip always convergent. 

In view of the failure of (0 for the second definition of 



CONCLUSION 

The extension of fluctuation theorems or the related 
non-equilibrium work relations to quantum systems re- 
quires the definition of work done on the system in a 
general non-equilibrium process. The aim of the present 
paper was to compare the implications of two possible 
definitions of work on the generalization of the Jarzyn- 
ski equation for a very simple, exactly solvable quantum 
system. This system is not in contact with a heat bath 
during the process and hence the work performed equals 
the change in its energy. 

If the system were described by classical mechanics the 
situation would be rather clear [3lj] . The energy differ- 
ence is given by the difference between the value of the 
final hamiltonian at the final phase space point and the 
value of the initial hamiltonian at the initial phase space 
point. Since the final state is a unique function of the 
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initial one the average in the Jarzynski equation is given 
by the sole equilibrium average over the initial condition 
of the system [l[ . 

If the system dynamics is described by quantum me- 
chanics at least two definitions for the work are conceiv- 
able. In the first one the the energy difference of the 
system is determined via measurements at both the be- 
ginning and the end of the process 3, Id, 15 , 27 1 . The 
energy difference then becomes a classical stochastic vari- 
able with transition probabilities derived from quantum 
mechanics and the corresponding average has to combine 
the average over the intial condition with the one over the 
transition probabilities. In this way the general validity 
of the Jarzynski equation can be established similar to 
the classical case described by a master equation [25| . 

If on the other hand the energy difference is described 
by the difference between the system Hamiltonian in 
Heisenberg representation at the beginning and at the 
end of the process it represents a genuine quantum ob- 
servable lESIi. In complete analogy with the 
classical case the average in the Jarzynski equation is 
then a sole average over the initial condition which is 
given by the trace with the initial density operator. This 
definition of work yields in the general case results dif- 
ferent from the first prescription. In particular we have 
shown that for a very simple, exactly solvable example 
the Jarzynski identity is violated if the second definition 
of work is used. 

The two definitions of work give identical results in the 
classical limit (3h — » as well as in the case where initial 
and final Hamiltonians commute with each other Tj^, 3] • 

At present it seems that the definition of work 
performed on a quantum system in a non-equilibrium 
process is to a certain extend a matter of taste. We 
therefore conclude that the question whether exact and 
general versions of non-equilibrium fluctuation theorems 
and of the Jarzynski equation applicable to quantum 
systems exist remains an open and challenging problem. 
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